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Abstract: Teak (Tectona grandis L.f.) is a popular hardwood species 
native to South and South-East Asia. The possible association of ampli¬ 
fied fragment length polymorphism (AFLP) marker with morphological 
variables of eleven important characters viz. girth at breast height (GBH), 
height of tree, bole height, branch knots, presence of fluting, spiral stem, 
leaf hair, leaf length/breadth ratio, branching pattern, bark colour and 
petiole shape for nine natural populations comprising 180 genotypes of 
teak (9 populations x 20 trees). The phenogram constructed using 
Euclidean distances for the eleven morphological characters showed that 
the populations were not grouped according to their geographical origin. 
The Mantel’s test for pairwise correlation between Euclidean distances of 
different morphological variables and genetic distances from AFLP data 
revealed that only petiole character (r =0.269; p =0.046) and height of 
tree (r =0.200; p =0.001) were significantly correlated with that of AFLP 
data matrix. The nine populations in this study covered a geographic area 
of about 1000 km stretch along the Western Ghat of South India. A test of 
correlation between genetic and geographic distance matrices revealed a 
significant positive correlation (r =0.475; p =0.009). The lack of perfect 
congruence between morphological and molecular data except for geo- 
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graphic distance, tree height and petiole character suggested that the 
morphological system might be useful for the morphotypes management 
but not appropriate to study the genetic structure of the teak populations. 

Keywords: AFLP, genetic, teak, morphology, Tectona grandis 

Introduction 

Teak has a wide but discontinuous distribution in India. It grows 
well from sea level to an elevation of 1200 m, and with an annual 
precipitation range from 800-2500 mm. Indian teak provenances 
show significant variation in tenns of adaptability, growth and 
form under a variety of edaphic and climatic factors. Gunaga and 
Vasudeva (2003) reported morphological variation in teak clones 
from Kerala and Karnataka states. Morphological variations in 
teak clones of Karnataka have been reported with respect to 
characters like floral traits (Vasudeva et al. 2001), vegetative and 
reproductive behaviors (Gunaga and Vasudeva 2003), seed ger¬ 
mination (Mathew and Vasudeva 2005), and early vigour 
(Mathew 2001). Bagchi (1999) measured five characters viz. 
total height, total bole height, diameter at breast height and 
crown width in 80 sets of teak trees selected from different plan¬ 
tations of Tamil Nadu, Kerala and Karnataka states. Kumar et al. 
(1997) considered height, diameter and basal area of different 
teak clones and reported significant genetic variation at the fam¬ 
ily level. Suri (1984) studied the quantitative and qualitative tree 
parameters in a provenance trial laid during 1931, comprising 
five seed sources from Burma, Kerala, Madhya Pradesh, Ma¬ 
harashtra and Karnataka. The statistical analysis showed that 
teak seeds from moist localities were bigger than those from dry 
zones. The tree height appears to be a function of site quality 
rather than origin of seed. Kjaer et al. (1996) while studying 
the genetic variation in teak from South Asia and Africa observed 
that differences between populations were much larger for mor¬ 
phological traits like stem fonn, relative branch size, wood den¬ 
sity and allozymes. 

Variability studies are important whenever selection forms the 
basis for genetic improvement (Bedell 1989). Kjaer et al. (2002) 
demonstrated the genetic and morphological variation in 
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Metroxylon sagu, a sago palm (a South -East Asian crop having 
high socio-economic importance) using AFLP markers. Sreeku- 
mar et al. (2006) conducted molecular and morphometric analy¬ 
sis of different Indian and Srilankan populations of two rattan 
species ( Calamus rivalis and C. metzianus) using RAPD markers 
and morphological characters. Multivariate analysis of RAPD 
and morphological data failed to separate rattan populations into 
two distinct species. Artocarpus altilis (Park) Fusberg is an im¬ 
portant fruit vegetable tree grown in homesteads of Southern part 
of India. In order to provide scientific management practices and 
conservation measures for this species, Sreekumar et al. (2007) 
investigated the pattern of morphological and genetic variation 
existing in six populations of this species using fifteen morpho¬ 
logical traits and AFLP markers. They found that genetic varia¬ 
tion from AFLP data is not reflected in quantitative morphologi¬ 
cal variables. However, the genetic and geographical distances 
were positively correlated which were well supported by the 
PCoA analysis and Dollo - parsimony tree. 

The main objective of teak plantation establishment in India is 
to produce high quality timber having acceptable growth rates. A 
large variation in growth and other economically important tree 
characters such as stem form, mode of branching, flowering 
habit and wood quality is observed in plantation-grown teak. 
Teak wood colour seems to be influenced by the site quality and 
moisture availability. As the golden brown wood colour is one of 
teak's most attractive timber qualities, the variation in wood col¬ 
our and texture has been widely studied (Kaosa-ard 1993). The 
success of planting programmes depends not only on the site 
quality but also on the genetic quality of the planting materials. 
Although the growth and yield of the plantations depend mostly 
on site selection and stem quality (i.e., straightness, persistence 
of stem axis and branching), is strongly controlled by genetic 
make-up (Keiding et al. 1986; Wellendorf and Kaosa-ard 1988; 
Kaosa-ard 1993). The results of international teak provenance 
trials carried out in Asian, African and Latin American countries 
have clearly demonstrated the influence of provenance or seed 
source on growth, stem quality, early flowering habit and health 
performance (Keiding et al. 1986). However, there is no avail¬ 
able information on relation between morphology and genetic 
data on teak. The objective of the present study was to investi¬ 
gate possible relations between morphology, genetic structure 
and geographical distribution of natural populations of teak from 
the Western Ghat region of South India and to identify possible 
AFLP marker associated with phenotypic characters. 

Materials and Methods 

Sample collection from natural forest 

The details of morphological characters from different study sites 
and leaf materials for DNA isolation were collected from nine 
natural teak growing forests of the Western Ghats belonging to 
the states of Kerala, Karnataka and Tamil Nadu (Table 1 & 2; Fig. 
1). These geographical areas were selected because they were the 
locations of origin of important provenances from the ‘very 
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moist’ to ‘moist’ teak forests of the Western Ghats identified for 
raising teak plantations. Historically, Konni, Parambikulam, 
Nilambur, Wayanad and Thrissur provenances were the sources 
of huge teak logs, and currently, majority of the teak plantations 
(70,000 ha approx.) and almost the entire seed production planta¬ 
tions of Kerala are located in the first four geographic areas 
(Prabhu, 2005). Teak from Nilambur Valley had been acclaimed 
as the most famous teak source used for shipbuilding and struc¬ 
tural purposes during 19 th century itself. Teak population sam¬ 
pled from Tamil Nadu state was from Indira Gandhi Wild Life 
Sanctuary located at the Eastern slope of the Western Ghats 
(Kala et al. 2005). Teak from Barchi (Haliyal Forest Division, 
Karnataka State) was famous as Dandeli teak because of the 
unique golden colour and fine grain (Dilip Kumar 2005). Virnoli, 
which was also part of Haliyal Forest Division was supposed to 
be the locality of origin of the famous Teli variety of teak which 
showed fast growth, shiny and smooth leaves that escape teak 
defoliator, Hyblaea puera attack (Katwal 2005). Moreover, ex¬ 
cept for Virnoli and Shimoga, these locations had very large, 
mature trees as well as smaller girth trees available for sampling. 
From each location, we collected expanding leaves and mor¬ 
phometric details from 20 randomly selected trees from a mini¬ 
mum of 2.5 km 2 area, sufficiently interior from reserve forest 
boundaries. 



Fig. 1: Map showing the location of natural populations of teak in the 
Western Ghats region selected for morphological and genetic analysis. 
WLS: Wildlife Sanctuary; IGWLS: Indira Gandhi Wildlife Sanctuary; 
P-V WLS: Peechi- Vazhani Wildlife Sanctuary 
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Table 1: Particulars of the site of teak populations sample collection from the Western Ghats used for AFLP analysis. 


Population 

Forest Division 

State 

Latitude (N) 

Longitude (E) 

Aerial distance from Konni (m) 

Konni (Konn) 

Konni Division 

Kerala 

09° 10' 

76° 57' 

- 

Thrissur ( Thri) 

Peechi- Vazhani Wild life Sanctuary 

Kerala 

10° 26’ 

76° 58' 

141 

Parambikulam (Para) 

Parambikulam Wildlife Sanctuary 

Kerala 

10°25’ 

76°45' 

141 

Nilambur (Nila) 

Nilambur South Division 

Kerala 

11° 18' 

76° 10' 

252 

Wayanad (Waya) 

Wayanad Wildlife Sanctuary 

Kerala 

11° 02' 

76°41' 

277 

Pollachi (Poll) 

Indira Gandhi Wildlife Sanctuary 

Tamil Nadu 

10 35' 

76° 52' 

158 

Shimoga (Shim) 

Shimoga Division 

Karnataka 

1381' 

75°58' 

539 

Barchi (Bare) 

Haliyal Division 

Karnataka 

15° 17’ 

74° 38' 

725 

Vimoli (Virn) 

Haliyal Division 

Karnataka 

15° 06’ 

74° 36' 

707 


Table 2: Characteristics of study sites of nine teak population from the 
Western Ghats. 


Population 

Mean 

Rainfall 

(mm) 

Tem¬ 

perature 

(°C) 

Alti¬ 

tude 

(m) 

Soil Type 

Konni (Konn) 

2229 

22-35 

220 

Sandy loam 

Thrissur (Thri) 

2838 

20-38 

55 

Oxisol/ red ferrallitic 

Parambikulam (Para) 

1890 

17-34 

560 

Sandy loam/laterite 

Nilambur (Nila) 

2666 

23-37 

68 

River alluvium 

Wayanad (Waya) 

2137 

17-34 

752 

Sandy loam/laterite 

Pollachi (Poll) 

2080 

14-34 

600 

Sandy loam/laterite 

Shimoga (Shim) 

1600 

20-35 

610 

clayey loam with alluvium 

Barchi (Bare) 

2700 

14-33 

573 

Deep sandy loam 

Vimoli (Virn) 

2700 

14-36 

599 

Deep sandy loam 


Morphological characters 

The following are the eleven morphological characters selected. 
Girth at breast height (GBH): the tree circumference measured at 
a point 1.37m above average soil level; height of tree: the verti¬ 
cal distance between ground level and tip of the tree; clear bole 
height: the distance between ground level and the crown point; 
fluting: the shape of the tree roots extending up from the butt of 
the tree, measured across the butt end of the log as the maximum 
inward deviation from roundness; spiral stem: growth of fibers in 
a spiral direction around the trunk of the tree causing twisting of 
timber; branching pattern: a secondary woody stem or limb 
growing from the trunk or main stem of a tree; knots on trunk: 
results from the inclusion of branches into the tree stem as the 
tree grows and increases in diameter; petiole shape: winged or 
sessile; leaf hair: hairs arise as the outward growth of epidermal 
cells; leaf ratio: the ratio between the leaf length and leaf breadth; 
bark colour: is the appearance of protective outer covering of the 
trunk. These characters were recorded in respect of each sampled 
tree. The details of the characters are listed in Table 3 according 
to the method suggested by Flush et al. (1972). 


Table 3: List of morphological characters used in the study 


SI. 

No. 

Character 

Measurements/Status 

1 

GBH 

in cm 

2 

Height of tree 

in cm 

3 

Bole height 

in cm 

4 

Branch Knots 

absence (0); Weak positive (1); Moderate positive 
(2); Strong positive (3) 

5 

Fluting 

absence (0); Weak positive (1); Moderate positive 
(2); Strong positive (3) 

6 

Spiral stem 

absence (0); Weak positive (1); Moderate positive 
(2); Strong positive (3) 

7 

Leaf hair 

absence (0); presence (1) 

8 

9 

Leaf length 

breadth ratio 

Branching pattern 

in cm 

Horizontal (0); Vertical (1); Horizontal to vertical (2) 

10 

Bark colour 

Greyish white (0); Brown (1); Greyish brown (2) 

11 

Petiole Shape 

clear (0); Winged (1) 


AFLP analysis 

Total DNA was extracted from 0.5 to 1 g of fresh leaf tissues 
using the modified CTAB protocol of (Doyle and Doyle, 1990). 
AFLP method was carried out following the standard procedure 
described by Vos et al. (1995). The reactions were carried out 
according to the manufacturer’s protocol (AFLP* Analysis Sys¬ 
tem 1 and Starter Primer Kit; Invitrogen Life Technologies, Inc., 
USA). The amount of DNA was adjusted to 250 ng for the TTcoRI 
and Msel digestion at 37 °C for 2 h. Ligation of TTcoRI and Msel 
adapters was done at 20 °C for 2 h. Five pi of the template DNA 
from a 1:10 diluted ligation mixture in TE (10 mM Tris-HCl, 1 
mM EDTA) buffer was used for pre-amplification with the 
pre-amplification primer mix. Twenty cycles of amplification 
were carried out at 94°C for 30 s denaturation, 56°C for 60 s 
annealing and 72°C for 60 S primer extension in PTC-100 Peltier 
Thermal Cycler (MJ Research, USA). Five pi of the 
pre-amplification product was used for selective amplification 
with ten primer combinations: E-AGG/M-CTT, E-AAC/M-CTG, 
E-AAC/M-CTT, E-ACT/M-CAQ E-ACC/M-CTT, 

E-ACA/M-CTT, E-ACC/M-CTA, E-AGG/M-CAG 

E-AGG/M-CTG and E-AGG/M-CAT selected out of 64 tested 
combinations. Primer labeling was perfolined by phosphorylat- 
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ing the 5' end of the EcoRI primer with y-p 32 ATP and T4 kinase. 
Mix 1 and Mix 2 components were used for selective AFLP am¬ 
plification in which. Mix 1 was composed of 5 pi of the labeled 
selected £coRI primer and 45 pi of Msel primer. Mix 2 was 
made according to the protocol supplied by Life Technologies 
with the AFLP kit. Selective amplification was carried out using 
5 pi of the pre-amplification product as the template DNA, 5 pi 
of Mix 1 and 10 pi of Mix 2. A touch down phase of thirteen 
amplification cycles was carried out initially, with first cycle at 
94°C for 30 s denaturation, 65°C for 30 S annealing and 72°C 
for 60 s primer extension. The annealing temperature was de¬ 
creased by 0.7°C each cycle during the following 12 cycles. The 
final 23 cycles were run at 94°C for 30 s, 56°C for 30 s and 72°C 
for 60 s. Samples were denatured by adding an equal volume of 
formamide buffer (98% formamide, 10 mM EDTA pFl 8, 0.05% 
bromophenol blue and 0.05% xylene cyanol) and heated at 94°C 
for 3 min. Five pi of each sample was run in a denaturing 6 per 
cent, 0.4 mm polyacrylamide gel (IX TBE buffer) in a sequenc¬ 
ing gel electrophoretic apparatus (Biorad, USA) at 40 A, 900 V 
for 2 h. The gel was then transferred onto Whatman No.l chro¬ 
matographic paper, dried at 80°C for 2 h on a gel drier (Biorad, 
USA) and exposed the gel to X-ray film (Kodak®) overnight. 

Data analysis 

Morphological analysis 

The straight forward method of qualitative two state characters 
recorded as + and - were scored as 1 and 0. Qualitative multistate 
characters recorded were scored as additive coding, with state (0) 
character undetectable, state (1) weak positive (+), state (2) 
moderate positive (++) and state 3 strong positive (+++) (Sneath 
and Sokal, 1973) (Table 3). 

Cluster analysis of morphometric data 
Cluster analysis was done for representing morphological dis¬ 
tances between populations using the NTSYSpc 2.1 package 
(Rolf, 2000). The morphological data matrices for eleven char¬ 
acters recorded from nine populations were transferred to 
STAND option to perfonn the standardization of data. The stan¬ 
dardized data were used to compute pair wise population dis¬ 
tance matrices based on Euclidean distance coefficient in 
SIMINT option and then subjected to unweighted pair group 
method with arithmetic average (UPGMA) clustering algorithm. 
Cluster analyses for 180 genotypes individually were also per¬ 
formed using EUCDIST (Euclidean distance) coefficient option 
in SIMINT module and then subjected to UPGMA clustering (r = 
0 . 86 ). 

PCoA analysis of morphological data 

Principal coordinate analysis (PCoA) was carried out for a better 
representation of morphometric distances between trees of all the 
populations using the option STAND, SIMINIT, DECENTER 
and EIGEN options using NTSYSpc version 2.1 (Rolf, 2000). 

AFLP data analysis 

AFLP data for teak trees from natural forests were utilized to 
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investigate the possible relation between genetic variation and 
morphological variation existing in teak populations AFLP bands 
were scored for presence (1) and absence (0) of bands. The data 
matrices were grouped into nine populations and analyzed using 
POPGENE, Version 1.32 package for dominant markers and a 
pair wise comparison of populations was made according to 
(Sreekanth et al. 2012). The Pair wise genetic distances (Nei 
1978) obtained were subjected for correlation studies (Table 4). 


Table 4: Nei’s unbiased measures of genetic identity (above diagonal) 
and genetic distance (below diagonal) coefficients 


Pop 

Konn 

Thri 

Para 

Poll 

Nila 

Waya 

Shim 

Bare 

Vim 

Konn 

*** 

0.928 

0.92 

0.919 

0.918 

0.926 

0.931 

0.916 

0.917 

Thri 

0.073 

*** 

0.937 

0.921 

0.924 

0.913 

0.913 

0.907 

0.902 

Para 

0.083 

0.064 

*** 

0.926 

0.902 

0.908 

0.901 

0.903 

0.9 

Poll 

0.083 

0.082 

0.076 

*** 

0.909 

0.906 

0.895 

0.889 

0.892 

Nila 

0.085 

0.078 

0.102 

0.095 

*** 

0.912 

0.904 

0.906 

0.908 

Waya 

0.076 

0.09 

0.096 

0.098 

0.091 

*** 

0.908 

0.909 

0.901 

Shim 

0.071 

0.09 

0.103 

0.111 

0.1 

0.096 

*** 

0.942 

0.943 

Bare 

0.087 

0.097 

0.101 

0.116 

0.098 

0.094 

0.059 

*** 

0.946 

Vim 

0.085 

0.102 

0.104 

0.113 

0.095 

0.103 

0.058 

0.055 

*** 


See Table 1 for population information 


Correlation of genetic and morphological distances using Man¬ 
tel s test 

Pair wise morphological distance matrix was created using 
EUCDIST (Euclidean distance) option in the SIMINT module of 
NTSYSpc 2.1 for each morphological character individually, as 
well as for the overall morphological data set. In order to find out 
possible relation between the morphological variation and ge¬ 
netic variation, product moment correlation coefficients were 
calculated between the genetic distance matrices and morpho¬ 
logical distance matrices for each character individually as well 
as for the overall morphological characters using Mantel test 
(Mantel, 1967). Statistical significance was evaluated with 1000 
permutations using TFPGA software, version 1.3 (Miller, 1997). 

Correlation between geographic data and genetic data 
The values of latitude and longitude recorded for each geo¬ 
graphic location of each of the nine natural teak populations were 
entered in a prescribed fonnat of GenAlex6 software package 
(Peakall and Smouse, 2006) for creating pair wise geographic 
distance matrix. The correlation test between the geographic 
distance matrices and the corresponding genetic distance matri¬ 
ces (as given in the above paragraph of AFLP data analysis) 
was performed using Mantel test (Mantel, 1967) and statistical 
significance was evaluated with 1000 permutations using 
TFPGA software, version 1.3 (Miller, 1997). 

Results 

Clustering and PCoA of morphological characters 

The phenogram constructed using Euclidean distances for eleven 
morphological characters showed that the nine populations were 
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not grouped according to their geographical proximity except 
Para and Thri (Fig. 2). There existed two main clusters, first 
main cluster formed by Konn, Nila, Shim, Waya, Bare, Vim and 
Tami and second main cluster comprised of Thri and Para which 
clustered separately. The cluster analysis of 180 trees from nine 
populations for similarity in morphological characters showed 
that the trees were generally scattered among different clusters 
without any geographical grouping based on their origin. 

The PCoA analysis of 180 trees based on eleven morphologi¬ 
cal characters failed to reveal any major separation in the pattern 
of morphological variations clearly indicating lack of distinct 
grouping of populations based on phenetic structure. None of the 
populations from neighbouring geographic regions tend to fall 
together in a given proportion of the two dimensional factor 
spaces. Some of the Para individuals were separated based on 
morphology (Fig. 3). 



Euclidean distances 


Fig. 2 : UPGMA dendrogram of morphological characters obtained using 
Euclidean distance between the nine teak population. 



Fig. 3: Principal coordinate analysis (PCoA) of 180 genotypes using eight morphological characters. 


AFLP analysis 

The standard genetic distance (D) unbiased for the sample size 
(Nei 1978) for all pair wise population comparisons were pro¬ 
vided in the Table 4 and the data were used to test the correlation 
between genetic data and morphological data. 

Mantel Test 

The mantel test did not reveal any correlation between Euclidean 
distance of overall morphological characters and AFLP genetic 


distance data (r =-0.1463; p =0.1740). The Mantel's test of pair 
wise correlation between Euclidean distance of individual mor¬ 
phological characters and genetic distance of AFLP data revealed 
that two morphological characters i.e. petiole character (r 
=0.2688; p =0.0460) and height of tree (r =0.2002; p =0.0010) 
were significantly correlated to AFLP data matrix. The test of 
correlation between genetic and geographic distance matrices 
also revealed a significant positive correlation (r =0.4748; p 
=0.0090) (Table 5). 
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Table 5: Mantel’s Correlation between distance value based on pair wise 
morphological features and the corresponding pair wise genetic distance 
value. 


SI. 

No 

First matrix 

(Morphological dis¬ 
tance) 

Second matrix 

(Genetic dis¬ 
tance) 

Mantels r 

Probabil¬ 

ity 

1 

Fluting 

AFLP data 

0.1566 

0.1650 

2 

Spiral 

AFLP data 

0.0000 

0.5350 

3 

Leaf hair 

AFLP data 

-0.0970 

0.7310 

4 

Leaf ratio 

AFLP data 

0.0355 

0.4850 

5 

Branching 

AFLP data 

0.0684 

0.3650 

6 

Bark color 

AFLP data 

-0.2101 

0.9030 

7 

Knots 

AFLP data 

-0.0970 

0.7230 

S 

Petiole 

AFLP data 

0.2688* 

0.0460 

9 

GBH 

AFLP data 

0.1003 

0.2650 

10 

Height of tree 

AFLP data 

0.2002* 

0.0010 

11 

Bole height 

AFLP data 

0.0803 

0.0460 

12 

Overall morphologi¬ 

AFLP data 

0.1463 

0.1740 


cal data 




13 

Geographic distance 

AFLP data 

0.4748* 

0.0090 


* Significant correlation (p <0.05) 

Discussion 

AFLP analysis was applied to reveal patterns of genetic variation 
in teak. The comparison between phenotypic and genotypic data 
depicted a low correlation. The nine populations used in this 
study covered a wide geographic range in latitudes and geo¬ 
graphical distances. The multivariate analyses of morphological 
data clearly indicated that populations from the neighboring lo¬ 
calities did not show any tendency to group together (Fig. 2 and 
3). The PCoA analysis supported the observation that the se¬ 
lected morphological characters scattered along 2D-scatter plot 
without a phenetic gap between each individual population. The 
detailed analysis of individual characters revealed that there is no 
clear distinction between nine populations based on quantitative 
and qualitative morphological characters analysed. A few of the 
Para individuals from Parambikulam Wildlife Sanctuary clus¬ 
tered together in one axis because of higher GBH which ranged 
from 2400 cm to 6520 cm. The discrimination among population 
variables might be strongly controlled by abiotic factors. In 
Baobab (Adansonia digitata) populations across different cli¬ 
matic zones of Benin (West Africa), Assogbadgo et al. (2006) 
observed significant correlation of morphometric characters with 
abiotic factors of the environment. 

The correlation between genetic distances and morphological 
characters revealed positive correlation of petiole shape and 
height of tree with AFLP data (Table 5). Mantel’s test of signifi¬ 
cance further gave an additional support and revealed that the 
two characters were significantly correlated with that of geno¬ 
typic data. Similar result was obtained in an AFLP marker study 
on genetic variation in a teak clonal seed orchard comprising 
1200 trees of 31 clones established in Southern Kerala in 1985 
(Sreekanth and Balasundaran, 2012). Mantel's test did not reveal 
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any correlation between overall morphology and genetic distance. 
The results of the Mantel’s test indicated that the genetic varia¬ 
tion detected by AFLP is not reflected in quantitative and quali¬ 
tative morphological characters measured. Ehara et al. (2002) 
performed RAPD analysis of 38 accessions of Metroxylon sagu 
from Malay Archipelago and found that none of the morpho¬ 
logical variables were useful in identification of specific geno¬ 
types. In the present study, the lack of correlation between ge¬ 
netic data and morphological variables was clearly indicated in 
the dendrogram and PCoA. 

One of the aims of the study was to find AFLP markers that 
are significantly more likely to co-occur with a morphological 
trait than expected by chance, that is, a marker that has a 
statistically significant association with the phenotypic trait. 
Ideally, we could not specifically reveal unique AFLP bands 
related to the phenotypic variables such as height of tree and 
petiole shape existed. Additional experiments such as QTL map¬ 
ping studies are needed to identify specific genes or genome 
regions that might have direct influence on the observed mor¬ 
phometric variation. 

The nine populations in this study covered a wide geographi¬ 
cal area. The most genetically similar populations were Vim and 
Bare and they were geographically very close by location, while 
the genetically distant populations of Bare of Karnataka State 
and Tami of Tamil Nadu State were far separated. The geo¬ 
graphical and genetic distances were positively correlated (r 
=0.4748; p =0.0090). Mantel's test of significance also proved a 
significant relation between geographic and genetic distance 
showing genetic divergence between distant populations. This 
significant correlation indicated that long distance dispersal of 
seeds had been rare. This was well supported by the UPGMA 
dendrogram and PCoA analysis (Sreekanth et al. 2012). 

The morphological clustering showed no geographical struc¬ 
ture. On the other hand, the molecular analysis showed that the 
genetic structure was related to the sample collection site. The 
lack of perfect congruence between morphological and molecu¬ 
lar data suggested that the morphological system might be useful 
for the morphotype management but was not appropriate to study 
the genetic structure of the teak population. The spatial structure 
of the genetic diversity may be related to the evolution of the 
species and the discordance between the morphological and mo¬ 
lecular structures might have resulted from similar selection 
pressures at different places leading to similar forms with a dif¬ 
ferent genetic background. 

Conclusions 

In a previous study of nine natural teak populations of the West¬ 
ern Ghats of South India strong association between genotype of 
natural teak populations and their respective geographic origin 
has been reported through AFLP marker study (Sreekanth et al. 
2012). About 75% of the total genetic diversity occurred within 
populations while only 25 % of the diversity was present among 
populations. The present study showed very little relation be¬ 
tween morphological variations and geographic distances be- 
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tween populations. Moreover, variation between populations in 
respect of most of the economically important phenotypical 
characters is not supported by genetic diversity revealed through 
AFLP markers. These findings show that ecological factors such 
as altitude, rainfall, soil characters, etc. also play a major role in 
phenotype of teak trees. Teak is a cross pollinated tree species; 
bole height, log size, tree form, wood colour, etc. are some of the 
economically important characters of teak trees. Hence, it can be 
concluded that for tree improvement of teak through phenotypic 
selection process, the number of plus trees selected from a popu¬ 
lation should be sufficiently large in order to exploit the large 
intra-population genetic variation. The results of these studies 
further emphasize the importance of progeny trial in teak tree 
improvement. 
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